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BACKGROUND OF THE INVENTION 

Field of the Invention : 

10 

[1] This invention relates in general to digital communication signal processing 
techniques, and more specifically to modulation and demodulation techniques in multiple 
sub-channel communications systems. 

15 Description of the Related Art : 

[2] Multiple sub-channel modulation and demodulation engines, and corresponding 
transmitters and receivers, which are implemented with digital signaling techniques are 
known in the art. One application of such a device is to communicate several sub- 

20 channels of information within a single, broader bandwidth, channel. Such systems 
employ orthogonal frequency division multiplexed ("OFDM") modulation, multi-carrier 
transmission ("MCT") modulation, and others modulation techniques. The general 
modulation approach in such systems is to consolidate the multiple sub-channels as sub- 
carriers in a composite signal in the base-band, and then frequency-shift the consolidated 

25 signal to the allocated carrier signal band, usually centered about a carrier frequency. 



Although, it is also understood that such systems can operate in the base-band without 
frequency shifting, and that such systems can be utilized in other electromagnetic bands, 
such as infrared and visible light. Multiple sub-channel modulation and demodulation 
devices employ various topologies and architectures. Among these are a particular 
5 variety that employ a combination of inverse Fourier transforms (or Fourier transforms) 
with polyphase filters and a software commutator. This is typically accomplished within 
a digital signal processor, however software implementations in other kinds of computing 
devices are certainly available. The basic design parameters of such devices, that are used 
for the purpose of modulation, are characterized by a number of input sub-channel 

1 0 signals, or sub-carrier signals, that are sampled at an input, or base-band, sampling rate, 
and, that are spaced apart in frequency by an input frequency spacing that is ultimately 
translated through to a composite signal. In addition, such devices are characterized by 
the composite output signal that is sampled according to an output, or composite, 
sampling rate. Respecting demodulation, the same characteristics apply, but are related to 

1 5 the input sampling rate of a composite signal, and the output sampling rate and frequency 
spacing of plural output sub-channel signals. These characterizations of are of vital 
importance when considering an overall system design and architecture. 

[3] The design of a polyphase filter for band limiting each sub-channel according to its 
20 pass-band and frequency spacing is critical to the design of an efficient data 
communications system. So too is the communications baud rate, which is directly 
affected by the filter design and sub-channel spacing, and which is important to the 
spectral response of the composite signal transmitted in the allocated channel 
Previously, those skilled in the art have understood that there were various constraints 
25 placed on the relationship between the input sampling rate, the sub-channel frequency 
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spacing, and the output sampling rate. These constraints have been considered 
problematic respecting optimum system design, especially where maximum performance 
in marginal signal conditions were sought. 



5 [4] More particularly, polyphase filters operate by multiplying selected phases, or 
samples, of a filter impulse response with samples of one or more the aforementioned 
input signals. Prior-art multiple channel polyphase filters have synchronized the selected 
phases of the filter impulse response with the positions of a commutator of the filter. In 
prior-art multiple channel polyphase filters, a given position of the commutator has 
1 0 corresponded uniquely to a predetermined phase of the filter impulse response. Those of 
ordinary skill in the art of polyphase filter design have even accepted that in multiple sub- 
channel polyphase modulators, the sub-channel sampling rate must be an integral multiple 
of the input sample rate, which also implies that the channel bandwidth, or frequency 
spacing, must be an integer multiple of the input sample rate. 

15 

[5] The tight limitations on the relationship between input sampling rate, frequency 
spacing, and output sampling rate were partially alleviated by the teachings of the present 
inventor, McCoy, in a prior US patent number 6,134,268 to McCoy for APPARATUS 
FOR PERFORMING A NON-INTEGER SAMPLING RATE CHANGE IN A 

20 MULTIPLE CHANNEL POLYPHASE FILTER (the £ 268 patent), the contents of 
which are hereby incorporated by reference thereto. The '268 patent teaches a multiple 
channel polyphase filter that includes a processing system for accepting and processing 
'M' input channels of data, each sampled at an input sampling rate, wherein 4 M' is a 
positive integer. The processing system is programmed to provide a commutator for the 

25 multiple channel polyphase filter, wherein the position of the commutator is decoupled 
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from the phase of a filter impulse response selected for the position, thereby allowing the 
multiple channel polyphase filter to be operated at a sampling rate that is a non-integer 
multiple of the input sampling rate. The processing system is further programmed to 
operate the multiple channel polyphase filter at the non-integer multiple of the input 
5 sampling rate to obtain a non-integer sampling rate change. Other embodiments and 
applications taught by the '268 patent include a multiple channel polyphase filter, a 
multiple channel modulation engine, a corresponding multiple channel demodulation 
engine, a multiple channel transmitter, and a corresponding multiple channel receiver. 
While all of these embodiments teach the decoupling of the input sampling rate from the 
10 output sampling rate, each still suffers from a limitation defining the sub-channel 
frequency spacing by a fixed relationship between the input sampling rate and output 
sampling rate. Thus, even in view of the teachings of McCoy in the '268 patent, 
communication system designers are faced with a constraints on sub-channel spacing that 
is often times unable to deliver optimum performance in practical applications. 

15 

[6] Thus, there is a need in the art for a modulation and demodulation engine 
applicable to multiple sub-channel systems that allows for the arbitrary specification of 
sub-channel frequency spacing with respect to input and output sample rates. 



20 



SUMMARY OF THE INVENTION 
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[7] The need in the art is addressed by the apparatus and methods of the present 
invention. The present invention encompasses multiple channel modulation bank that 
includes a programmable processing system for accepting plural channels of input data at 
an arbitrary sampling rate, or baud, rate, forming a composite output signal at an arbitrary 
5 sampling rate. Also, a multiple channel demodulation bank that includes a programmable 
processing system for accepting a composite input signal at an arbitrary sampling rate and 
yielding plural channels of output data at an arbitrary sampling rate. The sub-channels of 
data are spaced at an arbitrary spacing. While the prior art has restricted the relationship 
between the input sampling rate, output sampling rate, and the frequency spacing in some 
1 0 way, the present invention allows any rational relationship between all of these. 

[8] In particular, a modulator for receiving plural sub-channel signals that are sampled 
at a base-band sampling frequency and separated by a frequency spacing, and, for 
generating a composite signal, combining the plurality of sub-channel signals, that is 

1 5 sampled at a composite sampling frequency is taught. The modulator includes an inverse 
discrete Fourier transform coupled to receive the plural sub-channel signals and transform 
them into plural time domain signals. It also includes a multiple channel polyphase filter 
that receives the plural time domain signals and outputs a plurality of filter signals to a 
commutator. The commutator fractionally samples the filter signals at a rate defined by 

20 the ratio of the frequency spacing and a greatest common denominator between the 
composite sampling rate and the frequency spacing. 

[9] In a refinement of the foregoing modulator, the commutator fractional sampling 
ratio is derived as a ratio of integers thus allowing the relationship to be expressed as any 
25 rational number. In a further refinement, the modulator further adds a wireless modulator 
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that mixes the composite signal with a wireless carrier for use in wireless transmission. In 
a further refinement, the transform, the filter, and the commutator are implemented with 
executable software on a processor. In a further refinement, the processor is a digital 
signal processor. 

5 

[10] The present invention also teaches a modulator for receiving plural sub-channel 
signals that are sampled at a base-band sampling frequency and separated by a frequency 
spacing, and, for generating a composite signal, combining the plural sub-channel signals, 
that is sampled at a composite sampling frequency. The modulator includes an inverse 

10 discrete Fourier transform that transforms the plural sub-channel signals into plural time 
domain signals. The transform resolution is defined by the ratio of the composite 
sampling frequency and the greatest common divisor between the composite sampling 
frequency and the frequency spacing. Also, the transform decimation rate is defined by 
the ratio of the frequency spacing and a greatest common divisor between the composite 

1 5 sampling frequency and the frequency spacing. The modulator also includes a multiple 
channel polyphase filter that receives the plural time domain signals and outputs plural 
filter signals. The filter has an interpolation rate defined by a least common multiple 
between the base-band sampling frequency and the composite sample frequency divided 
by the base-band sampling frequency. Also, the filter decimation rate is defined as the 

20 filter interpolation rate times the base-band sample frequency and divided by the 
composite sampling frequency. Finally, the modulator includes a commutator that 
fractionally samples the plural filter signals at a rate defined by the ratio of the frequency 
spacing and a greatest common divisor between the composite sampling rate and the 
frequency spacing. 

25 
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[11] In a refinement of the foregoing modulator, the commutator fractional sampling 
ratio is defined by the decimation rate. In a further refinement, the modulator further adds 
a wireless modulator that mixes the composite signal with a wireless carrier for use in 
wireless transmission. In a further refinement, the transform, the filter, and the 
5 commutator are implemented with executable software on a processor. In a further 
refinement, the processor is a digital signal processor. 

[12] The present invention also teaches a demodulator that receives a composite signal 
that is a combination of plural sub-channel signals and that is sampled at a composite 

10 sampling frequency, and, that discriminates the plural sub-channel signals each at a base- 
band sampling frequency and separated by a frequency spacing. The demodulator 
includes a commutator that fractionally distributes the composite signal into plural filter 
input signals at a rate defined by the ratio of the frequency spacing and a greatest common 
denominator between the composite sampling rate and the frequency spacing. The 

1 5 demodulator also includes a multiple channel polyphase filter that has plural filter inputs 
to receive the plural filter input signals. The filter outputs plural filter signals to a 
discrete Fourier transform. The transformtransforms the plural filter signals, and outputs 
the plural sub-channel signals. 

20 [13] In a refinement to the foregoing demodulator, the commutator fractional sampling 
ratio is the decimation rate. In a further refinement, a wireless demodulator is added that 
is used to receive a wireless carrier signal, and to discriminate the composite signal 
therefrom to enable wireless reception. In a further refinement, the commutator, the filter, 
and the transform are implemented with executable software on a processor. In a further 

25 refinement, the processor is a digital signal processor. 
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[14] The present invention also teaches another demodulator that receives a composite 
signal that is a combination of plural sub-channel signals, and that is sampled at a 
composite sampling frequency, and, that discriminates the plural sub-channel signals, each 
5 sampled at a base-band sampling frequency and separated by a frequency spacing. The 
demodulator includes a commutator that fractionally distributes the composite signal to 
plural filter input signals at a rate defined by the ratio of the frequency spacing and a 
greatest common denominator between the composite sampling rate and the frequency 
spacing. The demodulator also includes a multiple channel polyphase filter that receives 

1 0 the plural filter input signals and outputs plural filter signals. The filter has a decimation 
rate defined by a least common multiple between the base-band sampling frequency and 
the composite sample frequency divided by the base-band sampling frequency, and the 
filter has a interpolation rate defined as the filter interpolation rate times the base-band 
sample frequency and divided by the composite sampling frequency. The demodulator 

1 5 also includes a discrete Fourier transform that receives and transform the plural filter 
signals, and outputs the plural of sub-channel signals. The transform resolution is defined 
by the ratio of the composite sampling frequency and the greatest common denominator 
between the composite sampling frequency and the frequency spacing. Also, the 
transform decimation rate is defined by the ratio of the frequency spacing and a greatest 

20 common denominator between the composite sampling frequency and the frequency 
spacing. 

[15] In a refinement to the foregoing demodulator, the commutator fractional sampling 
ratio is a ratio of integers that define any arbitrary rational number. In a further 
25 refinement, a wireless demodulator is added that is used to receive a wireless carrier signal, 
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and to discriminate the composite signal therefrom to enable wireless reception. In a 
further refinement, the commutator, the filter, and the transform are implemented with 
executable software on a processor. In a further refinement, the processor is a digital 
signal processor. 

5 

[16] The present invention also teaches a method of modulating plural sub-channel 
signals that are sampled at a base-band sampling frequency and separated by a frequency 
spacing, onto a composite signal that is sampled at a composite sampling frequency. This 
method includes the steps of converting the plural sub-channel signals into plural time 

1 0 domain signals by performing an inverse Fourier transform that has a resolution defined 
by the ratio of the composite sampling frequency and the greatest common divisor 
between the composite sampling frequency and the frequency spacing. Also, the 
transform has a decimation rate defined by the ratio of the frequency spacing and a 
greatest common divisor between the composite sampling frequency and the frequency 

1 5 spacing. The next step is filtering the plural time domain signals using a multiple channel 
polyphase filter to produce plural filter signals. The filter has an interpolation rate 
defined by a least common multiple between the base-band sampling frequency and the 
composite sample frequency divided by the base-band sampling frequency. Also, the 
filter decimation rate defined as the filter interpolation rate times the base-band sample 

20 frequency and divided by the composite sampling frequency. The last step is fractionally 
sampling the plural filter signals at a rate defined by the ratio of the frequency spacing and 
a greatest common divisor between the composite sampling rate and the frequency 
spacing. 
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[17] There is also taught a corresponding method of demodulating a composite signal 
that is sampled at a composite sampling frequency into a plurality of sub-channel signals 
each sampled at a base-band sampling frequency and separated by a frequency spacing. 
This method includes the steps of fractionally distributing the composite signal to plural 
5 filter input signals at a rate defined by the ratio of the frequency spacing and a greatest 
common denominator between the composite sampling rate and the frequency spacing. 
Next, filtering the plurality of filter input signals using a multiple channel polyphase filter 
to produce plural filter signals. The filter has an decimation rate defined by a least 
common multiple between the base-band sampling frequency and the composite sample 

10 frequency divided by the base-band sampling frequency, and the filter has interpolation 
rate defined as the filter interpolation rate times the base-band sample frequency and 
divided by the composite sampling frequency. The next step is transforming the plural 
filter signals, using a discrete Fourier transform, into the plural sub-channel signals. The 
transform resolution is defined by the ratio of the composite sampling frequency and the 

15 greatest common denominator between the composite sampling frequency and the 
frequency spacing. Also, the transform decimation rate is defined by the ratio of the 
frequency spacing and a greatest common denominator between the composite sampling 
frequency and the frequency spacing. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 
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Figure 1 is a functional block diagram of a wireless communications system 
according to an illustrative embodiment of the present invention. 

Figure 2 is a functional block diagram of a radio transmitter according to an 
illustrative embodiment of the present invention. 
5 Figure 3 is a functional block diagram of a radio receiver according to an illustrative 

embodiment of the present invention. 

Figure 4 is an architectural diagram of a prior art modulation engine. 

Figure 5 is an architectural diagram of a modulation engine according to an 
illustrative embodiment of the present invention. 
1 0 Figure 6 is an architectural diagram of a prior art demodulation engine. 

Figure 7 is an architectural diagram of a demodulation engine according to an 
illustrative embodiment of the present invention. 

15 

DESCRIPTION OF THE INVENTION 

[18] Illustrative embodiments and exemplary applications will now be described with 
20 reference to the accompanying drawings to disclose the advantageous teachings of the 
present invention. While the present invention is described herein with reference to 
illustrative embodiments for particular applications, it should be understood that the 
invention is not limited thereto. Those having ordinary skill in the art and access to the 
teachings provided herein will recognize additional modifications, applications, and 
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embodiments within the scope thereof and additional fields in which the present invention 
would be of significant utility. 

[19] An illustrative embodiment communications systems according to the present 
5 invention is illustrated in Figure 1. A plurality of terminal radio units 2 communicate 
through radio signals 4 with a radio repeater system 6. The radio terminal units comprise 
a transceiver circuit, controller, and user interface (not shown). The radio repeater system 
6 comprises a plurality of radio repeaters 8 that operate on discrete frequency allocations. 
Each radio repeater 6 comprises a transceiver. The transmitted and received radio signals 

10 4 are combined and coupled to one or more radiating structures with radio frequency 
combining circuits 10. The operation of the radio repeater system is generally controlled 
by controller 12, a computing device, the general function of which is understood by 
those of ordinary skill in the art. The controller 12 also comprises communications 
circuits. External networks, such as the public switched telephone network ("PSTN") 14 

15 and public or private data networks 16, are interfaced to the controller 12, which is 
operable to switch communications signals between the repeaters 8 and between the 
repeaters 8 and the PSTN 14 and networks 16. 

[20] More particularly, the illustrative embodiment in Figure 1 is implemented in a 
20 trunked land mobile radio system that employs FDM channelization of the allocated 
radio spectrum and TDMA access via packetized data for channel trunking management, 
system control, data communications, and voice communications within each of the FDM 
sub-channels. Such systems are deployed in the 800 MHz SMR band, the 900 MHz 
band, and certain other VHF and UHF bands. However, such techniques are equally 
25 applicable in any reasonable frequency band, as are those of the present invention. The 
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controller 12 may be any of a variety of computers, processors, microprocessor, or other 
suitable digital controlling devices as are understood by those of ordinary skill in the art. 
In Figure 1, three repeaters 8 are shown, however, those skilled in the art will appreciate 
that the number of repeaters is dependent upon the radio spectrum allocated to the 
5 system and may range from one to twenty, or more. The controller 12 provides various 
kinds of control of the resources within the repeater systems 6, including interconnecting 
radio and wireline communications resources, generating and interpreting communication 
protocols, encoding and decoding voice and data, and other control functions. 

10 [21] In an illustrative embodiment, the channel allocations are 25 kHz wide, and the 
aforementioned FDM channelization sub-divides each 25 kHz channel into four sub- 
channels, which are spaced 4.8 kHz apart. The TDMA channelization further divides 
each sub-band into two communications time slots with multiplexed data packets that are 
arranged in time to accomplish simplex, half-duplex, and duplex communications in real 

15 time. The illustrative embodiment system operates through the transmission and 
reception of four bit symbols at a 4 kHz symbol rate. The symbols are modulated and 
demodulated to the RF band according to a 16-QAM modulation scheme with a 16-point 
constellation in phase and amplitude, as is understood by those of ordinary skill in the 
art. The present invention teaches novel approaches to the base-band modulation and 

20 demodulation of such multiple sub-channelized systems. 

[22] Figure 2 is a functional block diagram of a radio transmitter 20 according to an 
illustrative embodiment of the present invention. A plurality of data signals 22, which 
may be arranged as a plurality of data channels, are coupled to a data/voice conditioning 
25 circuit 24. The data voice conditioning circuit may employ any of a variety of error 
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correction, interleaving, or other coding techniques that may be applicable to the 
communications objective of the systems. Such signal conditioning circuitry is 
understood by those having ordinary skill in the art. The voice data conditioning circuit 
24 is operative to generate 'M' base-band signals 25 corresponding to 4 M' sub-channels, 
5 that are sampled at a first base-band sampling rate. 4 M' is an integer value that is greater 
than one. In one illustrative embodiment, 'NT is equal to four. The data/voice 
conditioning circuit 24 is coupled to a processing system 26, which is a modulation 
engine, for accepting and processing the 'M' base-band signals. The processing system 
26 comprises a conventional digital signal processor ("DSP") 28 and a conventional 

10 memory 30 including RAM for providing an input buffer 34, and ROM for storing pre- 
programmed parameters and software, such as filter coefficients 32, a multiple channel 
polyphase filter program 36, and an inverse discrete Fourier transform ("IDFT") program 
38. Processed communication signals, in the form of frequency multiplexed composite 
signals, are output from the processing system 26 to a digital to analog converter 

1 5 ("DAC") 40, which converts the digital base-band signals to analog base-band signals. 
The analog signals output from the DAC 40 are coupled to an RF modulator 42, for 
modulating a wireless carrier with a frequency-multiplexed composite signal, which is an 
RF signal, having 'M' sub-channels. The signals are then radiated by antenna 44 for radio 
electromagnetic coupling via radio waves to one or more receivers. The DSP 28, memory 

20 30, and DAC 40 can be any of the variety of such devices known to be suitable for 
telecommunications applications by those or ordinary skill in the art, or that may later 
become known. Operation of the processing system 20 in accordance with the present 
invention will be described further herein after. 
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[23] An exemplary receiver 48 in an illustrative embodiment of the present invention 
appears as a functional block diagram in Figure 3. The receiver 48 is operable to 
demodulate a plurality of 4 M ? sub-channels that are received from the transmitter 20. 
The receiver 48 comprises an antenna 50 for intercepting radio signals. The antenna 50 is 
5 coupled to a conventional receiver front end RF demodulator 52 for receiving the radio 
signals and converting them to base-band analog signals. The RF demodulator 52 output 
produces a composite signal that is coupled to analog to digital converter ("ADC") 56. 
ADC 56 is operable to convert the analog signals into base-band digital signals, having the 
C M 5 sub-channels of information, into a composite signal having a sampling rate and 

10 derived from the £ M' sub-channels. The ADC 56 is coupled to a processing system 58 
for processing the signals and for controlling the receiver 48. The processing system 58 
comprises a conventional digital signal processor (DSP) 60 and a conventional memory 
62. The DSP 60, memory 62, and ADC 56 can be any of the variety of such devices 
known to be suitable for telecommunications applications by those or ordinary skill in the 

1 5 art. The memory 62 comprises software elements including executable code and variables 
storage for programming the processing system 58 in accordance with the present 
invention. The memory 62 includes filter coefficients 64 describing a filter impulse 
response in accordance with the present invention, as are understood by those of ordinary 
skill in the art. In addition, the memory 62 includes an input buffer 66 used for storing 

20 input data samples. The memory 62 further comprises a multiple channel polyphase 
filter program 68 and a discrete Fourier transform ("DFT") program 70, in accordance 
with the present invention. Operation of the processing system 58 in accordance with 
the present invention will be described further herein after. 
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[24] In several of the following illustrative examples, a multiple sub-carrier modulation 
system, developed by ComSpace Corporation and known in the marketplace by the 
"DCMA" trademark, is used to exemplify the prior art and the advantageous teaching of 
the present invention. The DCMA systems modulation approach transmits and receives 
5 a single or multiple sub-carriers within a FCC allocated radio channels that are spaced at 
25 kHz centers. Such channel allocations exist in the United States in the 800 MHz and 
900 MHz bands. However, the present invention is in no way limited by these examples, 
and, those of ordinary skill in the art will appreciate that the teachings herein are 
applicable to a vast array of communications systems and encompass a vast range of 

10 communications bandwidths. The exemplary DCMA systems modulates each sub- 
channel with a frequency offset that takes on one of four possible values with respect to 
the FCC allocated channel center frequency, and these values are; -7200 Hz, -2400 Hz, 
+2400 Hz, and +7200 Hz. Modulations is by 16-QAM comprised of 4-bit symbols 
transmitted at a symbol rate of 4000 baud. According to the aforementioned frequency 

1 5 offsets, the sub-carrier channel spacing is 4800 Hz. The output sampling rate of the 
combined sub-carriers ultimately transmitted over the allocated channel is 260 kilo 
samples per second ("ksps"). 

[25] Reference is directed to Figure 4, which is an architectural diagram of the prior art 
20 modulation engine as used in the DCMA modulation scheme. A plurality of ; M' sub- 
channel signals SM(m) 80 are input to upsample and pulse shaping filters 82. The filters 
82 are digitally implemented finite impulse response ("FIR") filters in one embodiment 
and output time domain signal PM(n) 84 to frequency shifting circuits 86 that shift the 
base-band signals according to the aforementioned frequency offsets. The frequency 
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shifted signals are summed be a summation function 88 and the combined output signal 
y(n) 90 is coupled to RF transmitting circuitry (not shown). 



10 
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[26] Mathematically, these operations can be expressed as follows. 

X 



0) 



[27] Where D is the polyphase filter decimation rate, lis the polyphase filter 
interpolation rate, m is the input signal time at the input signal sampling rate, and n is the 
output time at the upsampled output rate. And, s t (m) characterizes one symbol of input 
as a function of time. 

[28] Given that the upsample rate, establishes a modulo ratio, the expression becomes: 
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[29] And since any realizable digital filter must be truncated in time, the summation is 
limited to the number of taps per phase of the filter, the expression becomes: 
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[30] Adding the frequency shift and summing the plural singles, the output y(n) 90 
becomes: 
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[31] One of the advantages of the present invention is the reduction in processor 
operation required to achieve the desired modulation and demodulation of the plural sub- 
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channels. Respecting the prior art modulation engine of Figure 1, the number of processor 
multiplication operations (a reasonable indication of processor demand) required to 
implement the modulation engine is as follows. 



mult 




taps I phase 



+ Af(4 +/>))/„ 



(5) 



[32] Where p is the equivalent number of multiplies required in order to generate the 
complex sinusoid, and f so is the final sampling rate. Given the parameters of the DCMA 
system, the number of multiplies in the modulation engine of the repeater is: 



[33] Figure 5 is an architectural diagram of a modulation engine according to an 
illustrative embodiment of the present invention. The input signals SM(m) 100 are 

1 5 coupled to an inverse discrete Fourier transform ("IDFT") 102. The output bins of the 
IDFT are coupled to a multiple channel polyphase filter 104. The outputs of the filter 
104 are combined by a commutator 106 to yield the combined signal output 108. 
While the architecture in Figure 5 differs from the prior art architecture described with 
respect to Figure 4, the function of the modulation engine in Figure 5, in the general case, 

20 is mathematically the same. The structure in Figure 5 can be mathematically described, 
and the present invention further derived as follows. Consider a composite expression for 
the modulation engine, which combines Equations (1) and (4). 



Nmuit = (2 * 4 * 1 6 + 4 * (4 + /?))260000 - (37.44 + 1 Mp)MMPS 



(6) 




(7) 
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[34] Which is algebraically manipulated to yield Equation 8. 



x M 

y(n) = ^ h ( nD ~ m/ )2 s ii m y 2mihT * (8) 

m = -x l=\ 

5 [35] Use the substitution, n^r+sM, where r is an indicator of the commutator cycle 
position in time, and s is an indicator of overall time, to yield a modified summation. 



y(n) = y{r + sM) = ^ h(nD - ml)^ Sl (m)e jMr+sm > T > (9) 



10 [36] Equation 9 will be readily understood by those possessing ordinary skill in the art. 
An important advancement in the art made by the present invention occurs when the 
constraint that ^7; is equal to kiM, which equates^ ,to a ratio of integers, otherwise 
stated as a rational number. It is not required in general, but for illustration purposes, k is 
set to 1 and Equation (9) rewritten as: 



* M x M 

y(n)= ^ h(nD- mP^s^y 2 "^ 11 ™ = J h{nD-mt& s t ( m y , " t (10) 
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[37] Which is algebraically manipulated in the same fashion as between Equations (1) 
and (2) to yield: 



°° n {apsj phase -\ / 
y(n) = ^ h(nD - mf)s r (m) = ^ h(lm+ nDmod l)s r 
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[38] The commutator position s r {m) is interpreted as the rth IDFT output bin at time, 
m. In Equation (10), the sM term in the phasor repeats every 2jt radians and thus drops 
out of the equation. This filtering operation has the same number of processor multiplies 
as a single sub-channel rate change. The added processing cost for computing the multiple 
5 sub-carrier modulation bank is the cost of computing the IDFT, which performs the 

mixing operation. However, the IDFT is computed at the symbol rate and not the output 
sampling rate, and is therefore substantially less processor intensive. This demonstrates 
that the illustrative embodiment or the present invention not only achieves the desirable 
decoupling of the input frequency spacing from the filter parameters (decimation rate, 
1 0 interpolation rate, and number of phases) but, does so at a substantially reduced 
processor load. 

[39] Now, considering the forgoing advancement in view of the illustrative embodiment 
DCMA system, which operation parameters were enumerated herein before, it will be 

1 5 appreciated that the foregoing modulation model can be applied to the modulation bank 
for the downlink. The downlink is the communications path from the repeater systems 
to the terminal units. As the repeater systems must communicate simultaneously with a 
plurality of terminal units, it must do so on each of the sub-channels simultaneously as 
well. However, since the output sampling rate in DCMA is fixed at 260 ksps, it is 

20 necessary to contemplate the more complex sampling rate change this implies (4 ksps to 
260 ksps). Consider, generally, a process where the DFT is greatly over sampled in order 
to accommodate a difficult rate change. Such an operation can be mathematically 
constructed as follows. 
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[40] The first summation generally defining the filtering operation and the second 
generally defining the transformation. Again applying an advancement of the present 
invention, let^T^ be a rational number, a/b. Note that at this point in the design process, 
it is unknown how to select the IDFT resolution, M. Using this substitution, Equation 
(12) becomes: 



M 



y{n)= ^ h(nD- mfy^s^my 
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[41] Since^T; equals afb and this is k/M 9 then b equals aM/k, and k equals 1 , so b=aM, 
therefore Equation (13) becomes: 



M 

y(n)= ^ h{nD-miy} j s l {m)e ]2 "' alM 

/=1 



(14) 



[42] Again, applying the same practical application of a finite time filter, Equation (14) 
becomes: 



M ( 



y{n)= ^ h(mI+nDmodl)\s l 



m=0 



nD 

Tj 



■ m \e 



(15) 



[43] Where N tapsNyq is the number of taps per phase of the Nyquist filter, which is the 
filter chosen for use in the DCMA implementation of the present invention. The inner 
summation is the sub-sampled IDFT. A key difference between this equation and the 
previous equation is that the filter bank output is fractionally sampled. This is an 
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important advancement over the prior art. With this advancement, the baud rate can be 
related to the frequency by any rational number, and the output sampling frequency can 
be related to the frequency spacing by any rational number. This mathematical construct 
is very useful for any modulation engine or software radio. Some of the tradeoffs 
5 involved with such a structure do exist. The resolution of the IDFT is given by the 

output sampling rate divided by the greatest common denominator of the output sampling 
frequency and the sub-channel frequency spacing. The decimation rate on the output 
commutator is given by the sub-channel frequency spacing divided by the greatest 
common divisor of the output sampling frequency and the sub-channel frequency spacing. 

1 0 However, this advancement allows the designer of a communications system great 

flexibility in selecting the input sampling rate, the sub-channel frequency spacing, and the 
output sampling rate to suit system requirements. The variable factors which are 
adjusted to achieve this result are the number of terms of the IDFT, (the order of the 
IDFT), and the relationship between phase of the filter, (the sequence in which the 

1 5 phases are utilized). 

[44] Respecting the question of processor overhead burden of the DCMA illustrative 
embodiment of the present invention, in the worst case, the computational complexity of 
using this filter bank is as follows. 

20 

N mu it = 2 * N tapsl phase f so + C DFT f b (16) 

[45] Where C DFT is the number of processor multiplies required in order to compute 
the IDFT. 

25 
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[46] Again, for the DCMA illustrative embodiment, the following is an approximation 
of the complexity numbers. 

N mult = 2*16* 260000 + C DFT * 4000 = 8.32 + Q.0Q4C DFT MMPS (17) 

5 

[47] The question of how many processor multiplies will it take to implement the 
corresponding DFT. In this case there are far fewer inputs than total output points to be 
computed. Assuming the input is complex, the number of computations becomes: 

1 0 N mult = 8.32 + 0.004(4 * N sc * N DFT )MMPS (18) 

[48] It is important to note in this case that the number of operations required for the 
mixing operation can be increased above the required amount for the prior art 
implementation. This situation can arise if the number of sub-carriers is much smaller 
1 5 than the DFT size. However, in such situations, the utilized bandwidth in ratio to the 
output sampling rate will be small. If that is the case, the modulation bank can be 
implemented at a lower output sampling rate followed by an upsampler. The 
computational complexity of such an implementation in this case is as follows. 

20 N m «it =2*16 * 52000 + 4 * 4 * 65 * 4000 + N tap * 260000 - 5.82 + 0.26 * N tap (19) 

[49] Where N tap is the number of taps per phase to implement an up by 5 rate change. 
This second filter could be readily implemented in 40 taps per phase. This yields a final 
complexity of 16.22 MIPS, a savings of 20 MIPS over the prior art implementation. 
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[50] Respecting an implementation of the present invention modulation engine 
illustrated in Figure 5 to the DCMA illustrative embodiment, there are four input sub- 
channel sub-carriers SM(m) 100, each of which is sampled at 4 ksps. A total a 65 IDFT 
102 bins coupled to 65 Nyquist filter phases 104. The filter 104 interpolation rate is 65 
5 and the commutator 106 rate is 65 ksps. An upsample of the commutator output by 5 
brings the final output sample, y(n), 108 to 260 ksps. The filtering and commutation 
operation of this bank can be completely defined for all time in terms of three parameters, 
the filter phase, the data pointer position, and the commutator position. In the case of a 
DCMA illustrative embodiment 4 ksps baud system with frequency spacing of 4800 Hz 
1 0 and an output sampling rate of 52 ksps, these parameters are defined as follows. 



<p(n) - n mod 65 



(20) 



[5 1] The taps of the filter phase, are generated as follows. 



15 



\ =h(<j) + 65k), 



k = 0:N, 



tapsNyq 



-1 



(21) 



[52] The data pointer position is: 



20 



«(«)- 



n 



165 J 



(22) 



[53] The commutator position is given by: 



x(n) = (6n) mod(65) 



(23) 



24 



[54] The foregoing equation set does not specify the implementation of the IDFT. The 
IDFT has four inputs and 65 outputs. There are four active sub-carriers. The sub- 
carriers should are located at -7200 Hz, -2400 Hz, +2400 Hz and +7200 Hz. In order to 
5 accommodate these frequencies, the basic IDFT is modified as follows. 

jrM-J^m^ (24) 

[55] This additional shift is not expensive in terms of processor overhead because a 
1 0 direct implementation of the IDFT is used. The composite operation of the filter bank is 
summarized by the following steps: 

STEP 1: Compute the IDFT of one time-slice of symbols according to Equation (24). 
This will account for 65 separate filter histories, each at the symbol rate, (The original 
1 5 symbols will never be needed again.). 

STEP 2: Identify the filter phase according to Equation (20). 

STEP 3: Identify the data pointer position according to Equation (22). 

20 

STEP 4; Identify the commutator position according to Equation (23). 

STEP 5: Compute the dot product of the filter phase and the data vector according to the 
following equation: 
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tap^Nyq 

y(n)= gAW/i)+65*)-^ w («(n)-*) (25) 

[56] This operation is followed by a direct interpolation by 5 to achieve the design 
5 output sampling rate of 260 ksps, as discussed above. The interpolation by 5 operation 
is specified as follows. 



(j)(n) = n mod 5 

(j)(n) = n mod 5 (26) 



10 [57] The taps of the filter phase are generated as follows. 

\ = A(* + 5*), k = 0:N tapsNyq ~\ (27) 

[58] The data pointer position is: 



15 

fi(n)- 



[59] The output is computed as follows: 



(28) 



20 y{n)= 2M0(») + 5*)-;y(a(»)-*) (29) 



A.-0 
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[60] The teachings of the present invention are also applicable to the demodulation of 
multiple sub-carrier signals in a similar fashion to the previously discussed modulation 
aspects. For a more detailed understanding of the demodulation aspects of the teachings 
of the present invention, it is useful to consider the prior art, and especially the prior art 
5 DCMA systems. By way of review, DCMA modulation consists of a single (in the case 
of a terminal unit operating on a single sub-channel) or multiple sub-carriers with a 
frequency offset that takes on four possible values, -7200, -2400, +2400, or +7200 Hz 
from the allocated carrier frequency. The basic parameters of the system are a symbol 
rate is 4000 baud; the sub-carrier spacing is 4800 Hz, and an output sampling rate is 260 
10 ksps. 



[61] The prior art demodulation approach is diagrammed in Figure 6. A multiple sub- 
carrier signal Y(n) 1 10 is received from another circuit, which may be an analog to digital 
converter further coupled from an RF receiver demodulator. The signal is comprised of a 

15 digitally encoded signal at a relatively high sampling rate. A commutator 1 12 sequences 
the input signal through a number of positions equal to the number of sub-carriers 
represented in signal Y(n) 9 thus demultiplexing the received signals into their sub-carrier 
parts. Each sub-carrier part is frequency shifted by a suitable amount in Frequency Shift 
circuits 114. The frequency shifted circuits are output from the Frequency Shift circuits 

20 114 as base-band sub-carrier signals PM(n) 116. Each of these signals is sampled at the 
relatively high input 1 10 sampling rate. Down sampling to the sub-channel sampling rate 
and pulse shaping occur at Down Sampler/Pulse Shape circuits 1 1 8 for each sub carrier. 
The base-band signals, at the relatively lower output sampling rate are output at 120 as 
XM(rn). Note that 'tf is a representation of time at the relatively high input sampling 
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rate, and 'nC is a representation of time at the relatively lower output sampling rate. The 
total number of sub-carriers, or sub-channels is represented by the integer 'Af . 



[62] Mathematically, the processes illustrated in Figure 6 are expressed as follows. 

5 



[63] Where %i(m) is the output signal for each f sub-channel, D is the filter decimation 
rate, / is the filter interpolation rate, y(n) is the multiplexed input signal, f c is the carrier 
1 0 center frequency, Af is the sub-carrier offset frequency, and Ts is the sample time period. 
Then, using the substitution, n=r+sM, as was applied in novel aspects of the modulation 
engine calculations, and the summation for the novel aspects of the illustrative 
embodiment demodulation engine becomes: 




(31) 



x /( m )= 2 ^h{mD-(r +sAd)l)y(r+ sM)< 





(32) 



[64] Which is algebraically equivalent to: 



oc 



x i( m ) = 2 



e 



-j2m-(f c+ wyr s 




(33) 
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[65] In the flowing analysis, the inner summation, indicating the filter response, is 
separated out, and, the modulo count remainder of the decimation/interpolation process 
incorporated for thoroughness. 



5 y(m) = 2 h \ M1 



rnD-rl 
MI 



- 5 



+ (mD- rl)modMI \y{r + sM) (34) 
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[66] Using the substitution, s = 



mD-rl 



- s , and terminating the filter response time 



MI J 

according to the number of taps employed in the digital filter, the expression becomes: 



Ml 



y r {m) = ^ h(M£ + (mD-rl) mod M7);y r 



5 = 0 



mD-rl 



MI 



-s M\ (35) 
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[67] The operation of commutation is expressed by, y r (s) = y(r + sM), and is applied 
to Equation (35) to yield: 



MI 



yX m ) = 2 h{Mh +(mD- r/)mod Mi)y I 



~s=0 



mD-rl 



MI 



- s 



(36) 



[68] Equation 36 is the form of an ordinary polyphase filtering operation, with the 
slight modification of the rl term. This can be readily implemented in a digital signal 
20 processor, as will be appreciated by those of ordinary skill in the art. Finally, 



29 



substitution back in to Equation (33) yields the following complete expression of the 
output channels. 

M-\ M-\ 

x I {m) = y r ( m )= e J h s * s y r (™) (37) 

5 

[69] The discrete Fourier transform ("DFT") is not necessarily an ordinary DFT, 
although it can be in certain cases, depending on the relationship between Af and T s . In 
general, according to the advantageous teachings of the present invention, the relationship 
between Af and T s can be any rational number. In order to accommodate this possibility, 
1 0 the input to the DFT may have to be decimated. The processing structure employed to 
achieve the decimation, as well as the filter and transform according to an illustrative 
embodiment of the present invention is depicted in Figure 7. 

[70] Figure 7 is an architectural diagram of a demodulation engine according to an 
1 5 illustrative embodiment of the present invention. The composite multi-sub-carrier signal 
Y(n) 122 is input to a commutator 124, which commutes the signal 122 to the plurality of 
filter input signals 126. The filters 130 pulse shape the signals to YM(r,s,m) filter 
signals 132 which are input to discrete Fourier transform ("DFT") 134. The 'Af 
individual outputs xM(m) 136 are output from DFT 134, and are used for subsequent 
20 signal processing and utilization, as is understood by those skilled in the art. 

[71] Respecting the illustrative embodiment DCMA communications system, the 
input sampling rate is 260 ksps with 4 sub-carriers spaced at 4800Hz having a baud rate 
of 4ksps. The calculations yield values as follows; M=325,f c T s = -9/325, 4/7^=6/325, / 



30 



=1 and D=65. Consequently, the operation of the analysis bank can be expressed as 
follows. 



y T (s)-y{r+s32$ 



(38) 



325 



y r ( m )= J A(3255+(m65-r)mod32S)y r 



5=0 



m65-r 



325 



(39) 



and, 



324 



324 
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x t (m) = e Jl * trf325 ^e~ Jl2arU32S y (m) = ^ lS7r /325 ^^ 2jd ^ 6rmo<B2 ^ /325 c, 



y r (m) (40) 



r=0 



[72] Respecting the processor demand for this illustrative embodiment application of 
the present invention, the number of processor multiplications required in order to 
implement such an operation is approximated in the worst case as follows. 
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N mult = (2^ N taps + N dft )f h = (2 - N taps + 4N sc )f b 



(41) 



[73] The number of computations required for a 10 sub-carrier system assuming a 1 1 
symbol delay is 1 1 MIPS. The same system with a conventional implementation would 
20 require 70MIPS. On the receive side at least 3 taps of the filter must be sampled. A fine 
timing offset (early or late) may be realized by changing the phase of the filter. These 
additional outputs scale the complexity requirements. For 3 output taps, the complexity 
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requirements triple, which is true of both the prior art implementation and the illustrative 
embodiment implementation. 

[74] Thus, the present invention has been described herein with reference to a 
5 particular embodiment for a particular application. It is therefore intended by the 
appended claims to cover any and all such applications, modifications and embodiments 
within the scope of the present invention. 
Accordingly, 

10 WHAT IS CLAIMED IS: 
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